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NOMENCLATURE 

q  Uniformly  distributed  axial  load  per  unit  length 

P  ,  P^  Constant  end  loads 

L  Length  of  beam 

z  Axial  coordinate 

y        ;  Lateral  deflection 

t   ^\  ..  Time 

6   '  Variation 

T  Kinetic  energy 

V  Total  potential  energy 

m  Density 

E       •   Modulus  of  elasticity 

I  Moment  of  Inertia 

U  Strain  energy  of  bending 

W  Potential  energy  of  the  external  forces 

2 

CO  Separation  constant 

V 

X  =  —     ■  Non-dlmenslonal  axial  variable 

L 

P, 

B  =  -^  Ratio  of  end  load  P,  to  total  distributed  loads 
qL  1 

F  =  - —  Distributed  axial  loading  parameter 
EI 

n,  r  Subscripts  (Integers) 

f  Natural  frequency  of  vibrating  beam 

2  4 

„   CO  mL  ^.      , 

R  =  Eigenvalue 

EI 
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INTRODUCTION 

A  knowledge  of  the  effects  of  longitudinal  inertia  force  due  to  thrust 
on  the  vibration  and  stability  characteristics  is  important  in  the  design  of 
a  variety  of  aerodynamic  vehicles  and  their  components  especially  high  speed 
aircraft,  missiles  and  rockets,  where  inertia  and  friction  drag  forces  mani* 
fest  themselves  as  axial  loads.   In  recent  years  much  attention  has  been 
given  to  the  effects  of  axial  loads  created  by  the  thrust  on  the  natural 
frequencies  and  stability  characteristics  of  uniform  beams.   Timoshenko  [l] 
presented  the  natural  frequencies  and  buckling  loads  for  uniform  beams  sub- 
jected to  constant  end  loads.   Krieger  |  2]  and  Durgreen  [3]  determined  the 
effect  of  an  axial  force  on  the  free  vibration  of  hinged  bars  and  columns. 
They  considered  the  distance  constant  between  two  ends  instead  of  usual 
theoretical  assumption  that  the  load  on  the  beam  remains  constant.   Tu  and 
Handelman  (  4  |  considered  the  olgenvolue  problem  of  lateral  vibrations  of  a 
beam  under  initial  linear  axial  strees.   Real  \_5']   gave  some  results  of  the 
effect  of  longitudinal  acceleration  on  the  natural  frequencies  due  to  con- 
stant and  pulsating  thrusts.   Tyler  and  Rouleau  [6]  considered  an  analytical 
method  for  determining  the  deflection,  bending  moments  and  stabilities  of 
beam-columns  having  uniformly  distributed  axial  loads  by  the  Airy  integral 
function.   Oalley  ^7]  considered  the  natural  frequencies  of  vibration  of  a 
uniform  free-free  beam.   Recently,  Laird  and  Fauconneau  [8,  9]  have  investi- 
gated the  bounds  for  the  natural  frequencies  of  a  simply  supported  beam  and 
a  clamped  ends  beam  with  carrying  linearly  distributed  axial  loads.   The 
upper  bounds  have  been  obtained  by  the  Rayleigh-Ritz  method  and  the  lower 
bounds  by  Kato's  method  and  the  method  of  intermediate.   True  values  can  be 


predicted  from  the  upper  and  lower  bounds. 

In  this  report,  the  problem  of  stability  of  a  uniform  beam  with  linearly 
distributed  axial  loads  and  with  constant  end  loads  will  be  studied.   However, 
the  purpose  of  this  report  is  to  obtain  an  exact  solution  to  the  problem. 
Instead  of  finding  the  upper  and  lower  bounds.   The  natural  frequencies  for 
simply  supported  and  clamped,  as  well  as  for  the  free-free  beam,  are  obtained 
by  Frobenius*  method.   Finally,  the  results  are  compared  with  those  obtained 
by  Laird  and  Fauconneau  [9]. 


EQUATION  OF  MOTION  OF  BEAM  WITH 
UNIFORMLY  DISTRIBUTED  AXIAL  LOADS 


A  beam  of  imiform  cross  section,  subjected  to  constant  end  loads  P  ,  Pj 
and  a  uniformly  distributed  axial  loading  intensity  q,  is  investigated.   As 
shown  in  Fig.  1,  P  and  q  are  taken  to  be  positive  when  they  are  acting  in 
the  positive  direction  of  the  z-axis.   End  load  P-  acts  at  the  other  end  such 
that  the  force  system  is  in  equilibrium. 


Fig.  1.   Loading  Condition 

The  equation  of  motion  of  the  beam  can  be  derived  from  Hamilton's  prin- 
ciple, which  is  formally  expressed  by 


/ 

t 


(6T  -  6V)dt  = 


6/ 
t. 


(T  -  V)dt  =  0  , 


(1) 


where  T  represents  kinetic  energy,  V  the  potential  energy  of  the  system  and 
6  the  first  variation. 

Kinetic  energy  T  can  be  expressed  as 


T  =  1  /  m(-2)  dz  , 


(2) 


where  m  denotes  the  density  of  the  material  and  y  Is  the  lateral  deflection, 
of  the  beam.   V  can  be  broken  Into  two  parts:  the  potential  energy  of  the 
Internal  forces,  or  the  strain  energy  of  the  system,  expressed  as 

u  =  ^Ei  ;  (^)^dz  .  (3) 


2    h     ^z 


and  the  potential  energy  of  the  external  forces  W,  given  by 


W  =  -  i  /   (P  +qz)(^)^dz.  (4) 

2  Q      1         di 


Variations  of  T,    U  and  W  are  evaluted  as 


'^    '   ii'"r|f.^(6y)fa.  -l;\(|f)^z 


=   ji'"P|fa-t"^»^ 


=    /'"'lf.T»^»- 


1  L      .2  s2  2  1  L     N.2      2 

L    -2       .2 


z 


0    ^z^  ()z 


=     EI   /     T-^—^{bv)^z  , 


6W      =      -    /      1(P,    +qz)[t2  +^(6y)]^dz     +/     i(P,    +qz)(il)^dz 

q'-  x>  z      o  z  n  ^ 


and 


'i^'" 


6V  =  6U  +  6W 

L  .  2   V  2  L         V    \ 

'   =  EI  /  '^^fy(6y)dz  -/  (P,  +  qzi^ -^(6y)  dt  . 

Substitution  of  6T  and  6V  into  Eq,  (1)  yields 

t    L  2    2 

^H^   fm^^^Cay)  -  EI^f^(6y)  +  (P.  +qz)|^A(6y)]dz)dt  =  0.   (5) 


Integrating  by  parts,  the  first  term  of  Eq.  (5)  yields 


/V  ni^-^^(6y)dzdt  =  /  /  1  m^^^(6y)dtdz 
0  *^0 


L     ^v     t-    t,   -,2„ 

=  /   [n.<^(6y)]  ^  -  /  ^  in'i4(6y)dt  dz  . 

0     '■  t    t„   t    ^  t"^ 
0    0 


(6) 


The  first  term  of  Eq.  (6)  vanishes  automatically  because  variation  6y 
vanishes  at  t=t  and  t  ;  then 

*"!  ^   ?y  >  '^i  ^   A^y 

/V  ni<^^(6y)dzdt  =  -/  /  ni'^(6y)dzdt  .  (7) 

tO^'^^'^  tO^t 

0  0 

The  second  term  of  Eq.  (5)  yields  .  . 

t  L     2   ^2 

t  0    <^'  z  '^  z  .  ' 

0 


=  /'EijcH^^^y)^  -  /  ;Hf(6y)dz}dt 

'  ^''{f^rf  .■^<sy)]''  -  Crf^fiy)^''  -^  /''|^(6y)dz}dt .       (8) 

t^  l^^z  '^     0    ^z-*     0   0  ^z       ^ 


=  E 


=   A([(Pi +q2^-7<6y)]    -  /  ;7[(Pi +qz|7](6y)<i=sjdt  .  (9) 


and  the  last  term  of  Eq.  (5)  gives 

/  V  (P,  +qz)^|^(6y)dzdt 
0 

"0 
Substitution  of  Eqs.  (7),  (8),  (9)  into  Eq.  (5)  yields 

/  1/  In.  ^^  +  EI  ^  +  ^  [(P.  +  qz?-^])(6y)  dz  dt 
t  0  1  ^t^     ^^^   ^z  '  1     ^z'J 
0  (10) 

t       2         L       t       3       '  L 

+  /  1  [EI  ^1^  A(6y)]  dt  -  /  ^  [EI^(6y)  +  (P,  +  qzf^(6y)]  dt  =  0  . 

t      ^z^^z     0      t      ^z-*        ^     '^^     0 
0  0  .  . 

A  sufficient  condition  for  the  vanishing  of  the  left  hand  side  of  Eq.  (10) 
is  the  vanishing  of  each  part  of  the  corresponding  equation.   The  vanishing  of 
the  first  integration  for  arbitrary  6y  is  the  vanishing  of  its  integrand, 
which  yields  the  Euler  differential  equation: 

m  ^-4  +  EI  ^  +  ,—  I  (P,  +  qzl— J   =  0  .  (11) 


■■'S 


^j--^       ?^Z  ':'Z  '-   1       ?  Z-" 


The  vanishing  of  the  2nd  and  3rd  integrations  gives  the  required  natural 
boundary  conditions,  which  are  discussed  as  follow  :  r. 

(i)   Simply  supported  beam 

The  deflection  and  moment  must  vanish  at  the  end  points.  Therefore,  the 
boundary  conditions  are 

y  =  0 


■  « 


*     See  reference  [l5] 


■» 


(il)   Clamped  beam 

The  deflection  and  its  derivative  must  vanish  at  a  clamped  end.   The 
boundary  conditions  are   "  .   ' 

y  =  0 

(13) 

^  y 

(iii)  Free-free  beam 

The  moment  and  shear  must  vanish  at  a  free  end.   The  boundary  conditions 
become 

■'  ^  =  ° 

07. 
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SOLUTION  OF  EQUATION  OF  MOTION 

Equation  (11)  is  a  linear  partial  differential  equation  with  variable 
coefficients.   The  solution  is  assumed  to  be  of  the  form 

y(z,t)  =  X(z).T(t)  .  ■    ' 

Upon  substitution  into  Eq,  (11)  this  yields 

5 =  .  —5.  ^uj  (15) 

m         X  T    . 


The  expression  on  the  left  hand  side  is  a  function  of  z  only  and  the 

expression  on  the  right  is  a  function  of  t  only,  hence  they  are  both  equal 

2 
to  a  constant  co  .   Thus  Eq,  (15)  yields  two  ordinary  differential  equations, 

4  2 

— y  + |(P,  +  qz) —  I  -  =  0  (16) 

^2^   EI  dz  ^  1      dz-^    EI 

and 

d  T     2_ 

— r  +  uJ  T  =  0  (17) 

dt 

2 

The  type  of  transverse  motion  of  the  beam  depends  on  the  value  of  (J  . 

Thus  the  solution  to  Eq.  (11)  falls  into  the  following  three  cases: 

2  - 

(1)  to     >   0     The   solution  to  Eq.    (17)    is   T(t)   =  Acosujt   +  Bsln<x»t   . 

This    is   oscillatory  motion,   which   is   stable. 

2  _ 

(2)  UJ   =0  The  solution  to  Eq.  (17)  is  T(t)  =  At  +  B  .   This  gives 

the  static  buckling  case. 

2  - 

(3)  U!  <.  0  The  solution  to  Eq.  (17)  is  T(t)  =  Asinhtot  +  BcoshoJt  . 

This  represents  a  growth  of  T  when  t  increases,  which 
is  unstable. 
Introducing  the  following  non-dimensional  parameters 

z  1  qL  to  mL 

Eq,  (16)  now  becomes 

4 
d  X     dp,      dX-i 

__  +  F  — [  (B  +  X)— ]  -  Rx   =0.  (19) 

dx^     dx        dx 

In  view  of  the  difinition  of  parameter  B,  it  is  clear  that  for  a  given 
distributed  load  q,  the  following  cases  may  occur: 


9 


(1)  n  >  0  :   the  beam  is  entirely  in  compression. 

(2)  -1  <  B  <  0  :   the  beam  is  partly  in  tension  and  partly  in  compression* 

(3)  B  S  -1  :   the  beam  is  entirely  in  tension  since  the  tensile  end  load 
P  is  larger  than  the  distributed  loads. 

In  the  last  case,  the  problem  of  elastic  stability  does  not  exist. 

Determination  of  natural  frequencies  and  the  mode  shapes  involves  the 
solution  of  the  differential  equation  specified  by  Eq.  (19)  with  appropriate 
boundary  conditions.   By  the  well-known  method  of  Frobenius  [l2,13]  ,  a  power 
series 

r  2 

X  =  E  ax   =  a  +a,x+a,x  +.  .  .  (20a) 

r=0  ^       u    1     ^ 

can  be  assumed  as  the  general  solution  of  Eq.  (19),  Thus 

r-1  2       3 

X»   =  Z  a  rx     =  a   +  2a  x  +  3a  x  +  4a  x  +  •  •  •  ,    (20b) 
r=:0  '^  1^3      4 


0 

X"  =  E  a^r(r-l)x'""   =  [2]2a2  +  [^Va'^  +  WzV^  ""  '  ^^^""^ 


oo 

X*"  =  _S^  a^r(r.l)(r-2)x''-^  =  [sl^a^  +  [4]^a^x  +  [5]3a^x^  +  ..   (20d) 


r=0 


and 


CO 


Ifff 


,.0  r      ■   -■   -"'  =  W,a^-M,a3X.[6]^a^,2 


X   =  E  a^r(r-l)(r-2)(r-3)x     =  Wa%  +  C^la  x  +  [6]  a  x  +  •.   (20e) 


where  [r]^  =  r(r-l) (r-2).. . (r-n+1)  denotes  the  descending  factor. 
Substitution  of  Eqs.  (20)  into  lUq.  (19)  yields 


>V4  "-  f^^Vs''  ""  t^^Ve''^  ■"  ^^V?'*^  ■" 


■>,  . 


/  0  1 

+  F{a      +  2a   X   +  3a  x     +  4a,x     +  .    .    .  I        i^ 


^w  ^  ^ 


X   +  a^x     +  ax     +  •    • 


1=0. 


Collection  of   like  powers  of  x  results   in 

|[4]^a^  +FB[2]2a2   +  Fa^   -   Raj 


P\ 


+  ^ISJ^a^   +FB[3] 


2^   +  f[2] 


2^2   -^  2Fa2 


[6],a,    +  FB[4]-a     +  F[3],a     +  3Fa.    - 


J2  4 


2  3 


+  p\^7   +  FBL5]2a5   +  ^L^^^i,   +  ^Fa^ 


Ra 


Ra2|x 


Ra 


+  •    •    •  =     0   . 


10 


(21) 


Since  this  equation  must  hold   for  all  x,   the  coefficients  of   like   powers 

2 
of  X  vanish,  and    (n)   4n=n(n-l+l)=n  .     Thus 

[4]^a^  +  FB[2]2a2    +  Fa^    -   Ra^     =     0  ' 

[5]^a5   +  FB[3]2a3   +  F[2]^a2   -   Ra^      =     0         ■'■•'     Z;      ' 

■       .         [61    a.    +  FB[4],a,    +  F[3]^a,   -   Ra, 


■'4°6 


J2''4 


=     0 


[7]^a7   +  FB[5]2a5   +  f[4]   a^  -   Ra3      =     0 


[r]^a^  +  FB[r-2]2a^_2   +  ^C^^-^l   «r-3  '   ^^^.4     "     ° 


(22) 


11 


Substituting  the  coefficients  a  (k=4,5,6,»»»),  which  by  Eqs.  (22)  can 

be  written  in  terms  of  the  coefficients  a  ,  a  ,  a„  and  a  ,  into  Eq.  (20a), 

the  general  solution  of  Eq.  (15)  is 

3 
X  =  a^Xp  +  a^X^  +  82X2  +  a^X^  =  E  a^X^(F,B,R,x)  ,        (23) 

where  a  ,  a  ,  a.  and  a  are  arbitrary  constants  which  will  be  determined  by 
appropriate  boundary  conditions.  Also 

■•       R   4 
X  (F,B,R,x)   =  1+0+0+0+  -r^r-  x   +  •  •  •  ,  (24a) 

pV       Rv 

,.   X  (F,B,R,x)  =  X  +  0  +  0  -  — r-  + -^-  +  ...  .  (24b) 

2 

2       F'Br2]2  4   fF^I   5 
'X  X-(F,B,R,x)  =  X  +  0  -  — i=-^  X  -  -r4-x  +  .  .  .        (24c) 


and 


3   ,   fsC3]2  5   F(_-3]^  6 


X  (F,B,R,x)  =  X  +0 j^-=. —  X  -  ^  :^  X  +  •••  ,        (24d) 


The  recursion  relationship  for  Eqs.  (24)  is 

2  2 

for  n   4  . 


/ 
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SPECIAL  PROBLEM 

A.   Simply  Supported  Beam 

A  uniform  cross  section  beam  with  simply  supported  ends  is  considered* 
The  boundary  conditions  for  a  simply  supported  beam  are 

X(0)  =  X"(0)   =  X(l)   =  X"(l)  =  0.  (26) 

Substitution  of  Eq.  (23)  into  Eq.  (26)  gives 

a^  =  0  ,  (27a) 

"  '     a^  =  0  ,  (27b) 

a^X^(F,B,R,l)  +  a2X^(F,B,R,l)   =  0  (27c) 

and 

a^XJ(F,B,R,l)  +  a^X'^(F,B,R,l)   =  0.  (27d) 

In  order  to  obtain  a  nontrival  solution  of  Eqs.  (27),  the  determinant  of 
the  coefficients  of  a  ,  a  ,  a  and  a  must  vanish.   This  yields  the  characteris- 
tic equation 

X^(F,B,R,1).X^(F,B,R,1)  -  X^  (F,B,R,1).X'^'(F,B,R,1)   =  0  .      (28) 

For  given  values  of  F  and  B,  Eq.  (28)  is  a  polynomial  function  of  R  and 
was  solved  by  IBM  1620  computer  with  method  of  trial  and  error.   The  ratio  of 
convergence  was  quite  rapid;  no  more  than  100  terms  were  needed.   The  eigen- 
value R  in  Eq.  (28)  was  obtained  by  trial  and  error  for  given  loading  para- 
meters F  and  B. 

The  eigenvalues  for  the  first  five  modes  of  the  simply  supported  beam 
are  presented  in  Table  1.   The  Rayleigh-Ritz  upper  bounds,  the  lower  bounds 
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by  Kato's  method  and  the  method  of  Intermediate,  which  were  given  by  Laird 

and  Fauconneau  ^9],  are  also  Included  to  faclllate  comparison  of  these  methods* 

Analysis  of  Table  1  Indicates  that  the  Rayleigh-Rltz  upper  bound  is  the 
best  approximate  method.   For  B=0  and  1,  the  upper  bounds  actually  coincide 
with  the  results  obtained.   The  Rayleigh-Rltz  results  become  a  little  higher 
when  B  has  a  negative  value  and  F  has  a  large  value  and  the  beam  is  extremely 
close  to  buckling.   One  explanation  for  this  behavior  is  that  the  actual  mode 
of  a  simply  supported  beam  is  always  easy  to  assume  in  using  Rayleigh-Ritz 
method. 

The  effect  of  the  distributed  axial  loads  on  the  first  frequency  of  the 
simply  supported  beam  is  shown  in  Fig.  2.   The  ratio  of  the  first  frequency 
of  the  loaded  beam  to  that  of  the  unloaded  beam  is  taken  as  the  axis  of  the 
ordinate.   The  distributed  loading  parameter  F  is  taken  as  the  axis  of 
abscissa.   The  curves  are  plotted  with  various  values  of  B.   The  critical 
value  of  the  axial  distributed  loads  is  given  at  the  Intersection  of  the 
curves  with  the  horizontal  axis. 

B.   Clamped  Beam 

The  boundary  conditions  for  a  clamped  beam  in  non-dimensional  form  can 
be  expressed  as 

X(0)   =  X(l)   =  X'(0)   =  X«(l)  =  0  .  (29) 

Substitution  of  Eq.  (23)  into  Eq.  (29)  yields 

^       ^Q  =  0  1  .  (30a) 

a^   =  0  .  (30b) 

a2X2(F,B,R,l)  +  a^X^(F,B,R,l)   =  0  (30c) 
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and 

a^X^CF.B^R.l)  +  a^X'(F,B,R,l)  =  0.  (30d) 

From  the  set  of  Eqs,  (30),  the  characteristic  equation  for  clamped  beam 
is  derived  as 

X2(F,B,R,1).X«(F,B,R,1)  -  X^(F,B,R,1).XJ (F,B,R,1)  =  0  .     (31) 

For  given  values  F  and  B,  the  eigenvalues  R  are  calculated  by  trial  and 
error  as  in  previous  case. 

The  eigenvalues  for  the  first  four  modes  of  the  clamped  beam  are  listed 
in  Table  2.   The  results  obtained  by  the  Rayleigh-Ritz  method,  by  Kato*s 
method  and  by  the  method  of  intermediate  are  also  included  in  Table  2  to 
faciliate  the  comparison  of  these  methods. 

Examination  of  the  results  indicates  the  values  obtained  are  a  little 
less  than  the  Rayleigh-Ritz  method  but  are  much  larger  than  those  obtained 
by  other  methods.  Again,  the  Rayleigh-Ritz  method  gives  the  best  approximate, 

The  effect  of  the  distributed  axial  load  on  the  first  frequency  of  the 
clamped  beam  is  shown  in  Fig.  3,  which  is  plotted  in  the  ratio  of  the  first 
frequency  of  the  loaded  beam  to  that  of  the  unloaded  beam  vs.  loading  para« 
meter  F  with  several  values  of  B. 

C.   Free-free  Beam  ,     *^ 

In  this  section,  a  beam  of  uniform  cross  section  with  both  ends  free  Is 
considered.  The  boundary  conditions  for  a  free-free  beam  in  non-dimensional 
form  are 

k 

X"(0)  =  X"(l)  =  0  ,  (32a) 

X'"(0)  =  -  F(B  +  0)X»  (0)  (32b) 


24 


E 
01 

•o 

E 
«0 


<u 

u 

u 
o 

01 

3 

n 

01 
00 

u 


CM 

01 

n 
H 


o 
o 
ii 
n 


•o 

u 

c 

<0 

vt 

t-^ 

in 

o 

\0 

D 

•<-4 

E 

vO  <t 

00  in 

CM  r» 

-<  m 

»o  en 

o 

•D 

a> 

in  in  o  00 

CM  O  00  m 

00  00  \o  O 

o  oo 

00  o 

r>  — <  <o  <f 

X3 

t-i 

•  ■  •  • 

•  •  •  • 

•  •  •  • 

•   » 

•  • 

•  •  •  • 

^ 

o  (<-i  r^  M 

00  -^  in  o 

<f  \0  CM  vO 

o  o 

r^  00 

en  en  ^  O 

u 

o 

O  O  —1  <t 

en  r>-  — <  \0 

r^  m  <-i  o^ 

—  o 

O  — 

»j  >o  O  <t 
CN  00  \o  ^ 

2 

0) 

ui 

m  00  ^o  o 

«*  in  --  o 

en  m  \o  -< 

CO  — 

-<  en 

> 

jj 

a 

en  vt  o 

(^  <i  o 

d  n  00 

CO 

en  t^ 

CN  CM  vO 

q 

c 

^  CO 

'-*   fO 

^  m 

-<  CO 

-*  en 

»-) 

•^ 

1-1 

> 

C 
o 

CiO 

u 


Oj 


c   — 

O 


a 

3 


•o 

XI 

o 

C 

£ 

vt 

3 

iJ 

vo  <r 

o 

CD 

in  in  <o 

00 

J3 

•  •  • 

• 

O  en  r^ 

en 

^1 

O  O  -< 

vt 

o 

» 

in  00  \o 

O 

> 

O 

en  <t 

o> 

iS 

:2 

-^ 

en 

vt 

vO 

vt 

in 

in 

xO 

00 

• 

• 

• 

• 

o 

en 

r»- 

en 

o 

o 

f— < 

vt 

in 

00 

\f) 

o 

en 

<t 

O 

^^ 

en 

in  in  vD  00 

•  •  •  • 

O  en  1^  en 

O  C  —>  <t 

m  00  >o  o> 

en  <f  o 

«-■  en 


-^  CM  en  «3^ 


vt 

vO 

•-^ 

C> 

00 

ON 

en 

• 

• 

• 

• 

xO 

o 

o 

t»-i 

en 

IV. 

ft 

00 

vt 

in 

•— 1 

o 

en 

vt 

Ov 

en 

00 

rv 

<f 

en 

CM 

Ov 

• 

• 

• 

• 

00 

tN 

CM 

vt 

en 

rv 

CM 

00 

vt 

in 

t— 4 

O 

en 

vt 

O 

•— 4 

en 

in 

00  00 

<!•  en  CM 

c^ 

•  •  • 

• 

00  CM  CM 

^ 

cn  rv  CM 

•^  m  —• 

o 

en  <!■ 

o 

o 

o 

o 

•-4 

CM 

CM 

1^ 

• 

• 

• 

• 

vo 

o 

Ov 

Ov 

vO 

CM 

00 

00 

cn 

en 

in 

•— 1 

en 

en 

00 

•-I 

cn 

•-*  CM  en  »* 


CJv 

00 

vt 

t-i 

en 

C7> 

•-4 

• 

• 

• 

• 

in 

o 

vj- 

vt 

rv 

en 

CM 

CM 

en 

en 

VO 

tM 

en 

cn 

00 

t— 4 

en 

<M 

CJv 

vO 

•— « 

en 

Ov 

CM 

• 

• 

• 

• 

in 

Ov 

vt 

vt 

rv 

en 

CM 

CM 

en 

en 

vO 

CM 

en 

en 

00 

•-I 

CO 

vo  en 

--  rv  oo  in 

•  •  •  • 

o>  CM  rv  in 

oo  vD  ^  00 

CM  O  O  CM 

en  cn  rv 

—  m 


vO 

en 

00 

in 

en 

00 

<t 

• 

• 

• 

• 

o 

g 

in 

»— 1 

w-t 

CM 

vD 

en 

^^ 

f— 1 

en 

cn 

en 

rv 

r-t 

CO 

en 

vt 

CM 

in 

vt- 

ov 

rv 

• 

• 

• 

• 

o 

g 

in 

w—t 

t—4 

CM 

vO 

en 

^-1 

•— 1 

cn 

cn 

cn 

1^ 

^-1 

en 

CJV 

m 

vO 

m 

<t 

cn 

• 

• 

• 

• 

o 

vO 

CM 

:S 

o 

00 

rv 

CM 

rv 

vt 

en 

CM 

CM 

vO 

^•t 

O 

cn 

in 

en 

en 

<»■ 

O   OO 

• 

• 

•   • 

vt 

rv 

.t  vO 

vt 

vO 

CM  Ov 

CM 

oo 

vO  vj- 

CM 

CM  <0 

-H  en 

\D 

vo 

O 

en 

in 

■-< 

en 

• 

• 

• 

• 

:? 

rv 

in 

rv 

vO 

CM 

c^ 

CM 

oo 

vO 

vt 

CM 

CM 

vo 

■-) 

en 

«-•  CM  en  .J- 


o 
o 


•-<  CM  en  «;r 


o 
o 


•-*  CM  cn  <t 


o 
o 


o 

CM 


o 
en 


§ 


25 


XI 

0) 
4J 

• 

c 

<0 

in 

r^ 

<M 

lA 

r^ 

D 

•*H 

E 

^  o 

r-  00 

r^  CM 

^^  <ri 

—1  CO 

O 

•a 

o 

00  m  ^  o 

00  CO  00  CM 

\0  O*  \0  ITI 

o\  00  00  00 

m  «:t  i~»  o> 

J3 

•-I 

•  •  •  • 

•  •  •  • 

•  •  •  • 

•  •  •  • 

•  •  •  • 

J2 

m  m  <f  <M 

in  ^  O-  o 

fO  r«.  sf-  o> 

00  00  --  -- 

— <  <f  O  in 

U 

o 

r*  CM  c>  \o 

O  00  00  r- 

n  CO  f^  o 

•-<  00  1^  CM 

"^  <0  CVJ  CO 

01 

(U 

»H 

<-i  \o  O  m 

•->  (*1  ITI  \o 

-H  o  r>> 

O  O  \0 

o  o  m 

5 

i-) 

a 

<M  CM  in 

CM  — '  v* 

CM  «-"  fO 

CM  O  CO 

CM  O  CO 

9 

c 

•-<  m 

t-<  m 

—    CO 

•-•  cn 

-<  CO 

e 

a 


& 


0) 


o 

M-i 

0 

0) 
bO 


CM 


a 


c 

D 

o 
o 


o 


s. 


o 
o 
II 
(a 


<u 

D 

m 
> 

c 

60 


oi 


00  in 

CO  00  CM  CM 

•   •   •  • 

O  in  t-1  fsi 

•— •  O  0>  CM 

t-<  in  00  <f 

CM  •-•  in 

•-»  CO 


»D 

CO 

fO 

<}■ 

•* 

CM 

<f 

• 

• 

• 

• 

vO 

00 

CM 

o 

r~ 

CM 

CM 

CO 

f— 1 

vO 

^-i 

vO 

CM 

CM 

in 

•— t 

CO 

•O     JJ 

C       T^ 

in 

3    « 

in  <t 

0 

vt  in  vD  CM 

•  •  •  • 

\0  00  CM  — < 

r»  CM  CM  CO 

s.    ^ 

— "  vO  -J  \D 

CM  CM  in 

£•  >> 

— '  CO 

^   & 

o 

c-» 

^H 

w* 

*-^ 

c> 

<t- 

• 

• 

• 

• 

in 

vt 

o 

CO 

>o 

vt 

in 

CM 

CM 

^-1 

:? 

«-H 

CO 

o 

CM 

vO 

in 

CO 

ON 

CM 

• 

• 

• 

• 

VO 

p- 

t^ 

CM 

O 

co 

•—1 

lO 

f-^ 

CO 

vO 

(^ 

CM 

•— 1 

vt 

•—1 

CO 

r^ 

in 

CM 

in 

in 

•vJ- 

CO 

• 

• 

• 

• 

VO 

r^ 

00 

CO 

o 

00 

w-4 

vO 

•-4 

CO 

vO 

r~ 

CM 

t— 4 

vt 

t-1 

CO 

CO 

>* 

<r 

•-4 

• 

• 

• 

^-» 

CM 

c^ 

o 

o< 

in 

00 

in 

vt 

^-1 

O 

CO 

•-< 

CO 

m 

«-^ 

o 

CO 

CM 

00 

CM 

• 

• 

• 

• 

vt 

vt 

•— t 

(M 

CO 

vt 

•— 1 

O 

1—4 

wt 

00 

«M 

w— 

CO 

r-t 

CO 

vt 

in 

CO 

CO 

»t 

vO 

00 

• 

• 

• 

• 

vt 

vt 

CM 

CO 

CO 

vt 

f— 4 

a 

■—4 

^^ 

00 

CM 

1-4 

CO 

f* 

CO 

in  vt  i-i 

CO  vO  vO 

rv  .^  CM 

^  O  CO 
^  CO 


oo 


I  CM 

CO  r- 

•  •  •  • 

CM  O  tv. 

1  O  O  Ov 

f^  vt  o 

»H  o  CO 

"*  CO 


VO  CO 

in  CO  «;f  Ov 
•   •   •   • 

c^  in  o  rv 

t— (  CV  f— 4  <— * 

o  o  rv 

CM  •-"  CO 

'•*  CO 


CO 

V*  VO 
•-4  00 

•   • 

CM  O 

•-4  rv. 
O 

CM 


rv  00 

•  • 

Ov  O 

in  CO 

Ov  VO 

O  CO 

-<  CO 


fo 

o 

rv 

VO 

in 

CM 

rv 

• 

• 

• 

• 

o 

in 

1-1 

c^ 

4-4 

ov 

1-4 

i-i 

O 

o 

rv 

CM 

1-1 

CO 

1-1 

CO 

in 

CO 

•-4 

.-1 

1-4 

in 

m 

• 

• 

• 

• 

CM 

1—1 

o 

CM 

1-1 

rv 

VO 

CO 

o 

ov 

VO 

CM 

o 

CO 

•-< 

CO 

u 

"O 

u 
o 


i-l  CM  CO  «* 


CM  CO  «J 


1-1  CM  CO  >3' 


•-•  CM  CO  <t 


•-i  CM  CO  vS' 


o 
o 


o 
in 


O 
O 

o 


o 

o 


o 

o 

CM 


o 
o 


CO 

rv 


26 


C  CO 

3  -<  e 

O  T)     <U 

ti  •"  C   o 

o  c 


E 
« 


<u 

x: 


i-i 
o 


Q) 

3 

CO 


td 

<M 
V 

Ji 


o 

o 

II 

a 


3 
O 

a> 


■o 
o 

iJ 

0) 

,  E 


^     2 


3 
r-l 
(0 

> 

c 

bO 


ce: 


•n 

4J 

•«H 

3 

Cd 

o 

^5 

g. 

a 

>, 

CJ 

<§ 

O 

00 

>-i 

r^ 

-<  CM 

• 

• 

•   • 

vt 

0^ 

o  in 

m 

r-  «* 

o 

00  <}• 

CM 

O  m 

00 

r^ 

in 

.-< 

• 

• 

• 

o 

vO 

r~ 

vO 

r^ 

o 

»o 

fo 

a> 

<-^ 

O 

CM 

<-< 

fO 

vO 

ON 

1^ 

^0 

m 

O 

r- 

• 

• 

• 

• 

vt 

NO 

00 

in 

s 

o 

in 

CM 

o 

en 

^^ 

CO 

o 


CM 

vo  00  in 

•    •  • 

\0  On  in 

sj  O  <!■ 

O  On  in 

CM    O  CO 

— <  CO 


o 

•E 

O 


•-«    CM    CO    <t 


O 

o 


27 


0) 

u 

u 
o 


0) 
3 
rH 


C«4 


5 


11 
03 


T) 

iJ 

c 

«0 

CO 

00 

CM 

CM 

•-t 

3 

e 

vo  in 

O*  CM 

O 

r^ 

^0    «0 

o  r>. 

0 

■o 

01 

vo  r~  CM  r». 

tn  in  <M  sS' 

CM 

o  in 

P^ 

CM  CM  Vt  OO 

o  — •  m  <o 

ja 

>. 

UQ 

t-H 

•  •  •  • 

•  •  •  • 

• 

•  • 

• 

•   •   •   • 

•  •  •  • 

g 

J3 

r^  \o  <N  00 

fO  00  vO  <f 

00 

00  ON 

oo 

—  r-  CM  ^- 

CM  >0  CO  CO 

u 

p 

O 

O  in  r-  <f 

-<  o  cvj  in 

•-I 

in  r-- 

in 

CM  o  en  vO 

CO  NO  CO  <f 

01 

0) 

k 

^  «:r  00  vD 

rn  .-«  t-H  en 

CM 

r-»  CO 

o 

— •  «jf  NO  r~ 

CM  CO  CN 

> 

JJ 

CO  m  CO 

CO  c^  r^ 

CM  CM 

vO 

CM  — '  vt 

CM  —1  <T 

Q 

c 

^  m 

^  CO 

CO 

--   m 

-<  CO 

O 

o 


o 

4-1 

0) 

E 


O 

:2 


o 

3 
»-« 
(0 

^ 

to 

•t-4 


•H 

O 


bO 


U 
<U 

•a 
u 
o 


vt 

CJN 

00 

1^ 

CO 

<f 

<t 

• 

• 

• 

• 

s 

t— 1 

CO 

in 

in 

vO 

vt 

vt 

vt 

00 

NO 

CO 

CO 

00 

w-l 

CO 

C-l 

O 

NO 

rv 

CO 

CM 

00 

• 

• 

• 

• 

CM 

s 

rv 

in 

O 

00 

CM 

CO 

o 

o 

CO 

CO 

CO 

r^ 

•-• 

CO 

00 

r~~ 

CM 

o 

in 

CM 

• 

• 

• 

• 

CM 

o 

CM 

in 

c^ 

o 

On 

r-» 

w-i 

p~ 

(M 

c^ 

CM 

CM 

in 

t-l 

CO 

CO 

in  CM 

rv.  <f 

vJ- 

in 

•  • 

• 

• 

|v.  r» 

in 

vO 

O  in 

TV. 

in 

<r  >* 

00 

NO 

CO 

CO 

00 

—1 

CO 

c^ 

|v. 

NO 

fv 

r» 

<t 

NO 

• 

• 

• 

• 

CO 

On 

CM 

CO 

9—4 

O 

CO 

VO 

CO 

.— 1 

1-^ 

CO 

CO 

CO 

fv 

.-* 

CO 

in  r- 

<»■  >*  fv  O 

•    •    •    • 

00  O  00  o 

—<  NO  00  00 

CM  t*»  CO  O 

CM  CM  NO 

-'  CO 


Cl, 


o 
in  CM 

fv  <}•  <J-  \0 

•   •   •  • 

rv  rv  in  NO 

O  in  t~»  in 

.*  vt  00  nD 

CO  CO  00 

•-•  CO 


CO 

00 

oo 

rv 

rv 

in 

rv 

• 

• 

• 

• 

CO 

On 

CM 

00 

•-1 

O 

CO 

NO 

CO 

P~4 

1-1 

CO 

CO 

CO 

rv 

^>4 

CO 

o 

NO 

CM 

•J 

in 

ON 

CO 

• 

• 

• 

• 

00 

o 

CO 

o 

w~4 

NO 

00 

00 

CM 

rv 

CO 

o 

CM 

CM 

NO 

t-4 

CO 

rv  in 

m  CM  m  00 

•  •  •  • 

<t  m  CM  >-< 

fV  ON  00  .-< 

CM  <t  NO 

CM  •-<  vt 

—t  CO 


vt 

CM 

NO 

00 

o 

c^ 

rv 

• 

• 

• 

• 

CM 

m 

vt 

»-4 

CM 

CO 

in 

in 

•-4 

t-» 

o 

CM 

f^ 

vt 

w* 

fO 

CM 

m 

in 

in 

<t 

CO 

OO 

• 

• 

• 

• 

^-1 

ON 

vt 

o 

CM 

o 

vt 

On 

t-l 

vt 

NO 

rv 

CM 

•-< 

v* 

t-< 

CO 

c^ 

vt 

CM 

in 

Vt 

<7> 

• 

• 

• 

• 

CM 

00 

NO 

vt 

CO 

NO  vt 

rv 

CM 

CO 

CM 

CM 

f— 4 

vt 

f* 

CO 

fv 

vO 

CO 

in 

in  NO  «* 

• 

•  •  • 

t-> 

ON  vt  -< 

CM 

O  vt  Ov 

•— 1 

.d-  NO  r^ 

CM  —*  vt 

•-■  CO 

o 

CO 

CO 

NO 

oo 

NO 

• 

• 

• 

• 

CM 

00 

NO 

in 

00 

NO 

vt 

rv 

CM 

CO 

CM 

CM 

f-4 

J2i 

— <  CM  CO  vt 


-<  CM  CO  Vt 


»-<  CM  CI  «* 


— *  CM  CO  «;t 


^  CM  CO  vt 


o 


o 
o 


o 
in 


O 


CM 


28 


E 

n 

g. 
B 


u 
U 

o 

U-1 

0) 
3 
1-1 
TO 
> 
C 
Q) 
bO 


CM 

0) 

XI 
(0 
H 


pa 


•a 

XJ 

c 

<D 

D 

•rH 

E 

O 

•a 

<u 

X) 

J3     g 

t-H 

U 

XI 
O 

O 

01 

>H 

s 

4-> 

a 

o 

c 

.J 

•pi4 

■H 

D 
O 
XI 

u 
o 


•o 
o 

x: 

o 


(0 


3 
i-i 
CO 

(50 
•tH 

Id 


•a 
c 

O 
XI 


00 


s 


in 

«— t 

in 

in 

o 

^-4 

vO 

• 

• 

• 

• 

<N 

vt 

in 

in 

CM 

in 

00 

vfi 

o 

00 

vt 

CNJ 

o 

m 

—1 

en 

vO 

CO 

00 

c> 

• 

• 

• 

in 

<t 

vt 

VO 

r- 

O 

00 

in 

,_i 

«— 1 

O 

m 

*•*! 

m 

o 

,_( 

00 

vO 

<t 

o 

• 

• 

• 

• 

CM 

vO 

ON 

.— 1 

CM 

in 

O^ 

o 

o 

CO 

in 

CM 

o 

ro 

•—1 

n 

rsl 

cn 

CO 

r^ 

Ov 

o^ 

• 

* 

• 

• 

CM 

vn 

o 

«— 4 

fsl 

in 

C3N 

o 

o 

00 

in 

CM 

o 

ro 

•-" 

n 

CM 

O 

o 

in 

in 

r~l 

*— * 

• 

• 

• 

• 

o 

.-) 

in 

o 

•-< 

•—4 

Ov 

1— 1 

o 

r~- 

ro 

CM 

o 

CO 

t-4 

m 

•    vO 

00  fO  C3N 

I       •      •  . 

<t    CM  vO 

I  t-i  r>  CM 
00  v+  o 
•"-i  o  m 


o 

•^ 

vt 

00 

f— 1 

o 

CN 

• 

• 

• 

• 

o 

<f 

o 

vD 

•—4 

f— ( 

••t 

o 

00 

en 

CM 

o 

CO 

•-< 

en 

CM 

en 

vO 

00 

CM 

VO 

CM 

» 

• 

« 

• 

o 

vt 

O 

t^ 

*— 1 

«-4 

^H 

vi- 

O 

CO 

m 

CM 

o 

en 

•—t 

«n 

a> 

u 
o 


— '  CM  m  <j- 


— I  CM  m  vt 


u- 


in 

CM 


o 


in 

CM 


29 


E 
(0 

JO 


a> 
x; 


u 
o 


0) 

3 


CM 

o 
I 

II 


(0  CQ 


u 


CM 

0) 
t-l 

« 


■o 

<u 

c 

(0 

m 

in 

<t 

o> 

r^ 

D 

•.-4 

e 

ov  in 

O  S  on  O 

fo  —1 

CO  <N 

r^ 

in 

o 

•o 

d) 

—I  CO  --  r^ 

OO  »0  — '  vO 

\o  c^  cN  lA 

00 

r^  00  <r 

XI 

x>  g 

t-4 

•  •  •  • 

•   •  •  • 

•   •   •   • 

•  •  •  • 

• 

•  •  • 

x> 

<>  vo  n  00 

vO  r-~  i^  •-" 

0>  O  "-^  ^0 

CO  o  o  r^ 

CM 

yo  ->  <t 

CM  vt  <f 

u 

o 

vO  00  \o  o 

ro  vo  O  in 

»o  c^  o^  in 

c>  i~-  vo  in 

«N 

a> 

d) 

t^ 

v:t  >o  m  <r 

vj  in  •-"  o 

CO  ro  in  r-i 

(N  O  O  CM 

CM 

00  in  CO 

> 

u 

a 

o  <t  o* 

m  >*  <> 

m  fo  00 

CO  CO  r~ 

CM  CM  >0 

^ 

C 

'-•  r> 

— "  m 

-<  n 

—  CO 

!-•  fO 

n 
c 

3 
O 

o 


o 

JC 
0) 

o 


3 
t-l 
(0 
> 

c 

Q) 

CO 
*f-i 
M 


ce: 


(U 

•a 
o 


b. 


00 
vO  CM 

o>  \o  O  in 

•  •  •  • 

r>-  <t  <f  00 

lO  CO  \o  o 

<t  o  fo  in 

CO  <}■  o 

— '  CO 


o 

Ch 

CO 

vO 

in 

C7> 

• 

• 

• 

• 

o^ 

r^ 

o 

CO 

vO 

CO 

vO 

t— 4 

<t 

vO 

CO 

in 

CO 

vt 

a> 

•-I 

CO 

N 

•a   -u 

C        M-f 

<f 

3   oi 

0\  0^ 

o 

CO  \o  in  c> 

XI  ►.  x: 

•   •   •   • 

O  I^  O  CO 

VO  CO  \D  ^ 

K     ^ 

vj-  vo  CO  in 

ro  >;j-  o 

a   >, 

•-i  CO 

=    <s 

CM 

00 

<N 

\D 

CM 

00 

00 

• 

• 

• 

• 

^-) 

00 

p~ 

o 

CO 

in 

o 

VO 

Nt 

in 

o 

O 

CO 

vt 

O 

r-» 

CO 

vO 

vO 

O 

CO 

Vt 

f^ 

CM 

• 

• 

• 

• 

r» 

O 

c^ 

CM 

CO 

r~ 

•— * 

00 

vt 

in 

t— 4 

o 

CO 

vt 

o 

f-l 

CO 

rv 

VO 

t-^ 

fO 

«* 

rv 

CO 

• 

• 

« 

• 

r^ 

o 

o 

CM 

CO 

r~ 

f— 4 

00 

vt 

in 

•— 1 

o 

CO 

vt 

<3V 

W-* 

CO 

in 

f— 1 

00 

<!■ 

<f 

o 

• 

• 

• 

• 

|v. 

ro 

Ov 

Ov 

vt 

00 

CM 

CM 

CO 

CM 

in 

t-4 

CO 

CO 

00 

■-I 

CO 

r~ 

CO 

•—1 

CM 

O 

fv 

CO 

• 

• 

• 

• 

vt 

Ov 

vt 

Ov 

vt 

vt 

«-4 

vt 

CM 

c^ 

O 

—< 

CM 

CM 

rv 

f-l 

CO 

CO 

CM 

vO 

in 

CO 

o 

CO 

• 

• 

• 

• 

O 

•-I 

in 

CO 

rv 

CO 

•—1 

•-I 

CO 

CO 

VO 

CM 

CO 

CO 

CO 

•— 1 

CO 

in  CM 

in  CM  in  CM 

•  •  •  • 

ov  vD  CO  r^ 

o>  00  o  CO 

CM  O  -J  CO 

CO  CO  fv. 

•-I  CO 


CM 

f— 1 

in 

<r 

•— 1 

VO 

• 

• 

• 

• 

o 

t-H 

in 

CO 

r* 

CO 

t-^ 

^-1 

CO 

CO 

VO 

CM 

CO 

CO 

00 

— 1 

CO 

rv 

vO 

CO 

in 

CO 

00 

ov 

• 

• 

• 

• 

Ov 

VO 

fO 

rv 

Ov 

00 

o 

CO 

CM 

o 

t— 1 

CO 

CO 

CO 

rv 

1-^ 

CO 

»-<  CM  CO  <f 


•-»  CM  CO  «t 


»-<  CM  CO  Vt 


<-l  CM  CO  >* 


O 

o 


o 


o 


o 


in 

CM 

CO 

o 

in 

CM 

r- 

• 

• 

• 

• 

VO 

in 

vO 

—1 

—1 

<N 

in 

00 

f-i 

VO 

CM 

o 

CM 

CM 

vO 

<— 1 

CO 

CO 

^^ 

CO 

CO 

• 

• 

• 

• 

CO 

vt 

in 

CO 

(M 

CO 

00 

in 

CM 

00 

in 

vt 

CM 

CM 

VO 

ft 

CO 

CO 

00  o  00  -^ 

•  •  •  • 

CO  in  in  in 

CM  CO  CO  in 

CM  00  in  <j- 

CM  CM  VO 

<-"  CO 


CM  CO  <J 


o 

00 


30 


B 
n) 

0) 
.Q 


g. 


^1 

O 
(u 

0) 

n 
<u 

•fH 


in 

CM 


« 


c 

(0 

3 

•^ 

E 

O 

•o 

o 

U3 

0) 

.-1 
0 

GJ 

0) 

)-l 

> 

u 

a 

0 

c 

J 

•^ 

T) 

o 

c 

j:: 

3 

OJ 

o 

0) 

^ 

>> 

a 

u 

£i 

M 

(U 

» 

> 

0 

^ 

TO 

<D 
3 
t-i 
TO 
> 

c 
to 


Oi 


o 


00  f-t 


iTl  iTl    O 

lO  iTl    O 

CM  in  o 

CM  —<   s* 


vt 

fO 

vO 

f— < 

cn 

<!■ 

• 

• 

• 

• 

CM 

r»- 

o 

ro 

<t 

r~ 

vD 

vo 

<— 1 

in 

O 

in 

CN 

CM 

in 

•-4 

CO 

fn 

in 
in 

O 

Ov 

t— 1 

00 

CO 

• 

• 

• 

• 

<r 

m 

CO 

in 

in 

f— 1 

CM 

vO 

c^ 

in 

vO 

CVJ 

•-H 

<t 

.-1 

c^ 

t 


in 

00 

o 

<N 

r^ 

r-» 

in 

• 

• 

• 

• 

00 

CO 

o 

vt 

r^ 

vO 

t— 4 

1— < 

CM 

CM 

C-J 

1-1 

vt 

•— 1 

CO 

r^ 

3     « 

o^  in 

o 

vO  <f  CM  <r 

•^   &^•= 
».-^i? 

•  •  •  • 

OJ  r^  •— '  in 

«d-  r^  vo  o 

g.     X 

—'  in  O  in 

CM  CM  in 

c.   > 

—I  m 

--  CM    M   <]• 


-H  CM  m  «d- 


o 
o 


o 

CM 


CM  cn  <!• 


o 


31 


ft*       cr 

•^^    ^ 
"^^ 

II 

II 

CQ 

Cl, 

B 

(U 
CO 


-^    s. 


<M 
O 

>> 

u 
c: 


<0 

c 

(U 

e 
a 
•o 


0) 


o 
n 
"S 


o 

iJ 

O 

0) 


00 


o 

0=J 


NO 

o 


o 


o 


o 


oi:;eH  Aouanbaaj 


32 


and 


X"'(l)  =  -  F(B  +  l)X'(l)  .  (32c) 


(i)  For  the  general  case  when  B  t^  0   and  F  ^  0,  substitution  of  Eq.  (23) 
into  Eqs.  (32)  yields 

a^  =  0  ,  (33a) 

agX'^(F,B,R,l)  +ajX'^(F,B,R,l)  +  a^X^  (F,B,R,1)  +  a^X|J(F,B,R,l)  =  0,  (33c) 

ao^o(F,B,R,l)  +  a^Xj'(F,B,R,l)  +  a2X2(F,B,R,l)  +  a^X^"(F,B,R,l) 

=  -  F(B  +  l)[aQXQ(F,B,R,l)  +  a^X^ (F,B,R,1)  +  a2X^ (F,B,R,1) 
+  a2X«(F,B,R,l)]  .  (33d) 

Substitution  of  Eqs.  (33a)  and  (33b)  into  Eqs.  (33c)  and  (33d)  gives 


2i 

"0"0  '  ''3'-"3    FB 


a.X'J  +  a.rX'J  -  — ^  X"]   =  0,  (33e) 


and 


■  "fT  ^^^  "^  ^^^1^   =  °  •  <33f) 

Thus   the  characteristic  equation  for  free-free  beam  is 

M3. 

X"  X     +  F(B  +  1)X; i|   X,    +  F(B  +  1)X'  1 

^  jpB^l  1-" 

r  W3         ^_ 

"  ^^^3   "  li~  ^1  ^L  "^0  ""  ^^^  •"  ^^^0^       =     °   •  (34) 
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(ii)  When  B  =  0,  the  boundary  conditions  of  Eqs.  (32)  become 

X"(0)  =  X"(l)  =  0  ,  ^  (35a) 

X  (0)  =  0  (35b) 

and  * 

X  (1)  =  -  FX«(1)  .  (35d) 

Substitution  of  Eq.  (23)  into  Eqs.  (35)  yields 

32  =  0  ,  (36a) 

a^  =  0  ,  (36b) 

aQX'^(F,B,R,l)  +  ajXl^>(F,B,R,l)  =  0  (36c) 

and 

aQ[x|J'(F,B,R,l)  +  FX^(F,B,R,1)  +  a^[xJ|'(F,B,R,l)  +  FXJ(F,B,R,1)]  =  0.   (36d) 

The  characteristic  equation  is 

X'^[X"J' +  FXJ]  -  X'^[X|;' +  nc»^]  =  0.  (37) 

(iii)  IThen  F  =  0,  B  =  0,  the  boundary  conditions,  Eqs.  (32),  become 
'  ^   '       X'^(O)  =  X"(l)  =  X'"(0)  =  X"'(l)  =  0.  (38) 

Substitution  of  Eq.  (23)  into  Eq,  (38)  gives 

^2  =  0  ,  (39a) 

33  =  0  ,  (39b) 

a^XJj(F,B,R,l)  +  a^X'^(F,B,R,l)   =  0  (39c) 

and 
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aQX'^'(F,D,R,l)  +  a^x'"(F,B,R,l)   =  0.  (39d) 

The  characteristic  equation  is 

X"X'"  -  X"X"'  =  0  .  (40) 

0  1    10 

The  third  case  in  this  report  is  to  present  eigenvalues  for  a  free-free 
beam  by  the  same  method  as  those  obtained  for  a  simply  supported  beam  and  a 
clamped  beam.   Translational  and  rotational  modes  corresponding  to  a  zero 
natural  frequency  will  exist  if  the  beam  is  free  on  both  ends.   Using  a 
strict  definition  of  stability,  such  motion  would  be  unstable.   However,  the 
type  of  instability  in  which  failure  occurs  due  to  structural  damage  caused 
by  bending  is  of  primary  importance.   For  this  reason  rigid  body  translation 
and  rotation  will  not  be  considered  to  be  unstable  in  a  free-free  beam  case. 

The  results  are  presented  In  Table  3,   The  effect  of  the  distributed 
axial  loads  on  the  first  frequency  of  the  free-free  beam  is  again  shown  in 
Fig.  4,  as  before.  When  F  =  0  and  B  =  0,  this  case  reduces  to  a  uniform 
free-free  beam  without  axial  load.  Meirovltch  [l2]  gave  the  same  answers, 
R  =  500.564  and  3803.54,  for  the  first  and  second  modes. 
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Table  3.   Eigenvalues  of  the  free-free  beam 


B  =  0 

Eigenvalue  R 

F 

1st  mode 

2nd  mode 

3rd  mode 

0.0 

500.564 

3803.54 

14617.6 

5.0 

375.001 

3530.10 

14149.7 

10. 0 

249.721 

3255.28 

13680.8 

15.0 

138.039 

2980.76 

13211.3 

20.0 

57.2011 

2708.73 

12742.0 

22.5 

28.5197 

2574.50 

12507.7 

25.0 

5.27823 

2441.95 

12273.6 

B  =  +0.5 

Eigenvalue  R 

F 

1st  mode 

2nd  mode 

3rd  mode 

1.0 

450.910 

3694.51 

14430.7 

3.0 

350.721 

3475.86 

14056.5 

5.0 

250.080 

3256.63 

13682.0 

7.0 

151.181 

3037.09 

13307.3 

9.0 

59.5318 

2817.60 

12932.5 

10.0 

19.4391 

2708.02 

12745.1 

B  =  1.0 

Eigenvalue  R 

F 

1st  mode 

2nd  mode 

3rd  mode 

1.0 

426.031 

3639.95 

14337.2 

3.0 

275.375 

3311.70 

13775.9 

5.0 

124.763 

2982.47 

13214.0 

6.0 

52.2076 

2817.70 

12933.0 

6.5 

17.9028 

2735.32 

12792.5 

6.75 

1.53088 

2694.15 

12722.2 
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Tabic  3.   Eigenvalues  of  the  free-free  beam 


D  =  -0. 

25 

Eigenvalue  R 

F 

1st  mode 

2nd  mode 

3rd  mode 

1.0 

488.122 

3776.26 

14570.9 

10.0 

371.672 

3528.15 

14148.2 

30.0 

160.316 

2975.04 

13202.2 

50.0 

71.5717 

2462.97 

12261.3 

70.0 

14.0831 

2033.98 

11345.3 

73.0 

5.74743 

1977.44 

11211.2 

B  =  -0.5 

Eigenvalue  R 

F 

1st  mode 

2nd  mode 

3rd  mode 

1.0 

500.498 

3803.50 

14617.6 

20.0 

476.067 

3790.63 

14608.6 

40.0 

421.064 

3753.82 

14582.1 

60.0 

368.555 

3698.09 

14538.7 

80.0 

329.750 

3629.76 

14480.6 

■^v- 
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SUMMARY 

In  this  report,  a  direct  method  of  obtaining  the  natural  frequencies  of 
a  uniform  beam  carrying  linearly  distributed  axial  loads  is  proposed.  The 
governing  differential  equation  with  varying  coefficients  has  been  solved 
directly  by  Frobenius'  method.   The  solution  is  in  series  form.   It  is  not 
necessary  to  approach  these  problem  by  an  indirect  approximate  method  in 
seeking  the  upper  and  lower  bounds.  The  eigenvalues  of  the  simply  supported, 
clamped  and  free-free  beam  carrying  linearly  distributed  axial  loads  have 
been  presented  in  Tables  1,  2,  3,  for  a  range  of  loading  parameters  F  and  B* 

It  should  be  mentioned  that  the  equations  derived  in  this  report  are 
based  on  the  small  deflection  theory  by  which  it  can  be  considered  that  the 
axial  loads  remain  constant  in  magnitude  during  the  beam's  motion  and  that 
the  supports  are  free  to  slide  in  the  axial  direction.   The  equations  are 
also  valid  only  for  the  lower  modes  because  the  effects  of  shearing  deforma« 
tion  and  rotatory  inertia  become  [l3,14]  as  the  mode  number  increases. 

The  method  presented  in  this  report  could  be  extended  to  the  considered 
tion  of  beams  with  other  boundary  conditions  and  to  other  problems,  giving 
rise  to  differential  equations  with  variable  coefficients,  such  as  the 
determination  of  natural  frequencies  and  buckling  loads  of  beams  of  varying 
cross  sections.   Information  of  this  nature  would  be  valuable  to  designers  of 

high  speed  aerospace  vehicles. 

"  *        •■  .        ,"%■ 

.!'■'.      1  ■■'■<■'  ''  '  '         '%  -^ 
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In  this  report  the  stability  of  a  flexible  uniform  beam  with  linearly 
distributed  nxinl  loads  and  with  constant  end  loads  is  considered,   Thl8 
problem  is  reduced  to  that  of  finding  the  eigenvalues  of  a  fourth  order 
ordinary  differential  equation  with  variable  coefficients.   By  Frobenius* 
method  the  natural  frequencies  of  vibration  of  the  beam  for  bending  modes 
are  found  as  functions  of  the  loading  parameters.   The  numerical  values  of 
natural  frequencies  are  evaluted  by  trial  and  error.   The  fundamental  fre- 
quency vs.  loading  parameter  curves  Is  plotted  for  the  simply  supported, 
clamped  and  free-free  beams.   Values  of  the  critical  distributed  axial  loads 
are  found  at  the  intersections  of  these  curves  with  the  horizontal  axis 
where  the  natural  frequency  equals  zero. 


